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MIT  Lincoln  Laboratory  is^l4ilfi|ating  the  detection  and  identifica¬ 
tion  of  stationary  ground  targets  in  high  resolution,  fully  polarimetric, 
synthetic  aperture  radar  (SAR)  imagery.  This  article  (t)  provides  a 
brief  description  of  the  Lincoln  Laboratory  SAR,  (2)  describes  an 
optimum  polarimetric  processing  technique  used  to  construct 
minimum-speckle  SAR  intensity  imagery,  and  (3)  presents  examples 
of  polarimetrically  processed  imagery.  ®  1993  John  Wiley  S  Sons. 
Inc. 

I.  INTRODUCTION 

The  Lincoln  Laboratory  SAR  is  a  fully  polarimetric,  TT-GHz 
synthetic  aperture  radar  sensor  (1|.  The  polarimetric  capabili¬ 
ty  of  the  radar  is  used  to  enhance  the  quality  of  the  imagery 
taken  from  a  small  aircraft;  the  synthetic  aperture  permits 
data  to  be  processed  to  a  resolution  of  1  ft  x  1  ft  at  a  slant 
range  of  7  km.  The  sensor  was  developed  to  provide  a  high- 
quality  database  of  clutter  and  target-in-clu  ter  imagery,  for 
use  in  evaluating  the  performance  of  station  iry-target  detec¬ 
tion.  discrimination,  and  identification  algorithms. 

An  example  of  the  quality  of  imagery  gathe'ed  by  the  SAR 
is  presented  in  Fig.  1.  This  synthetic  aperture  radar  image  of  a 
golf  course  located  near  Stockbridge,  NY,  has  undergone  the 
optimum  speckle-reduction  processing  described  here.  Figure 
2  shows  a  close-up  photograph  of  the  golf  course;  note  the 
pond,  flag  pole,  putting  green,  and  line  of  four  trees  which  are 
visible  in  both  Figs.  1  and  2.  Because  of  its  high  resolution, 
the  Lincoln  Laboratory  SAR  can  resolve  individual  trees  and 
bushes,  as  well  as  the  pond  and  the  putting  green  shown  in  the 
image.  Note  the  1  ft  x  1  ft  resolution  permits  one  to  discern 
very  small  objects  such  as  the  flag  pole  located  in  the  center  of 
the  putting  green.  This  image  was  obtained  under  clear 
weather  conditions.  However,  the  quality  and  resolution  of 
the  SAR  image  would  not  be  degraded  in  the  presence  of 
dense  fog  or  thick  cloud  cover.  Thus  a  SAR  sensor  has  a 
significant  advantage  over  optical  sensors;  the  image  quality  is 
not  dependent  on  weather  conditions,  and  the  sensor  can  be 
used  either  during  the  day  or  at  night. 

Figure  3  shows  the  Lincoln  Laboratory  SAR  Sensor  and 
lists  some  of  the  system  parameters.  The  SAR  Sensor  is  an 
airborne,  instrumentation-quality  radar  carried  on  a  Gulf- 
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stream  G1  aircraft.  The  radar  antenna  is  housed  in  a  specially 
designed  radome  mounted  beneath  the  aircraft.  The  radar 
transmitter,  dual-channel  receiver,  and  digital  recording 
equipment  are  all  carried  by  ihe  airplane. 

Since  a  principal  objective  of  the  Lincoln  Laboratory 
studies  is  to  evaluate  the  benefits  of  fully  polarimetric  radar 
data  for  stationary-target  detection  and  identification,  polari¬ 
zation  purity  is  essential.  To  achieve  polarization  purity,  a 
special  corrugated  horn  antenna  with  a  Fresnel  lens  was 
designed,  providing  very  pure  horizontally  and  vertically 
polarized  transmit  waves.  The  radome  was  designed  to  mini¬ 
mize  cross  coupling  between  the  horizontal  and  vertical  po¬ 
larizations.  The  radar  transmits  horizontal  and  vertical  polari¬ 
zations  on  alternate  pulses;  dual  receiver  channels  measure 
both  retii.-ns  simultaneously.  Inertial  velocity  estimates  com¬ 
pensate  foi  aircraft  motion  between  the  horizontal  and  verti¬ 
cal  transmit  pulses. 

An  in-scene  corner  reflector  calibration  array  comprised  of 
several  high-quality  trihedrals — and  dihedrals  oriented  at  0, 
22.5.  and  4.5" — is  used  for  polarimetrically  calibrating  the 
imagery.  The  polarimetric  calibration  scheme  is  described  in 
ref.  2.  During  each  flight,  data  are  gathered  and  digitized  in 
real  time  with  a  28-channcl  Ampex  recorder.  The  data  are 
then  brought  to  the  Lincoln  Laboratory  ground  processing 
facility,  where  S.AR  image  formation  is  performed.  Special- 
purpose.  high-speed  digital  processing  hardware  is  used  to 
construct  the  imagery  and  perform  the  polarimetric  cali¬ 
bration. 

A  SAR  is  a  radar  that  synthesizes  a  long  aperture  as  an 
aircraft  flies  along  its  path.  Thus  a  SAR  can  achieve  cross¬ 
range  resolutions  that  could  otherwise  be  attained  only  with  a 
long  antenna.  In  SAR  mode  the  Lincoln  Laboratory  radar  has 
1  ft  X  I  ft  resolution.  Range  resolution  is  achieved  by  using 
6(K)-MHz  bandwidth  pulses.  To  achieve  1  ft  azimuth  res¬ 
olution,  a  synthetic  aperture  of  approximately  150  m  length  is 
constructed  by  processing  1  s  of  data  as  the  plane  flies. 

SAR  processing  can  produce  high-resolution  images,  but 
the  process  is  subject  to  a  considerable  amount  of  speckle  in 
the  images  because  of  the  cohe'ent  nature  of  the  imaging 
process.  Noncoherent  spatial  averaging  of  high-resolution 
pixel  intensities  can  be  used  to  reduce  image  speckle.  For 
example,  we  have  significantly  reduced  image  speckle  by 
averaging  4x4  pixel  clusters  of  single-polarimetric-channel 
intensity  data  into  effective  1-m  x  1-m  resolution  pixels.  How¬ 
ever,  the  speckle  reduction  was  obtained  at  the  cost  of 
degraded  image  resolution.  A  new  technique  described  here, 
the  polarimetric  whitening  filter  (PWF),  uses  a  polarimetric 
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Figure  1.  SAR  Image  of  a  golf  course  in  Stockbridge,  NY  (1  ft  x  1  ft  resolution).  The  sensor  was  flown  at  an  altitude  of  2  Km  with  a  look-down 
(depression)  angle  of  22  5°.  giving  a  slant  range  of  7  Km.  PWF  processing  was  used  to  produce  this  minimum  speckle  image.  The  radar  is 
located  at  the  top  of  the  image  looking  down,  therefore,  the  radar  shadows  go  toward  the  bottom  of  the  page 


methr.'d  of  speckle  reduction  that  preserves  image  resolution 
(.T4|.  This  filter  processes  the  complex  (HH.  HV.  VV)  data 
into  full-resolution  pixel  intensities  in  a  way  that  minimizes 
SAR  image  speckle.  This  method  is  based  on  a  mathematical 


model  that  characterizes  fully  polarimetric  radar  returns  from 
clutter.  By  using  this  polarimetric  clutter  model,  we  can 
derive  an  algorithm  that  shows  how  fully  polarimetric  data  can 
be  combined  into  minimum-speckle  imagery. 


Figure  2.  Optical  photograph  of  one  side  of  the  golf  course.  This  photograph  was  taken  on  November  t99t.  Some  visible  features  in  the 
photograph  are  the  flag  pole  in  the  center  of  the  putting  green,  the  open  area  pond,  and  the  four  large  trees  adjacent  to  the  pond. 
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Figure  3.  The  advanced  detection  technology  sensor  (right).  The  sensor  platform  is  a  Gulfstream  G1  aircraft  shown  in  flight  (left).  The  ADTS 
radome.  located  at  the  bottom  of  the  aircraft,  was  built  by  the  Loral  Corporation.  System  features;  SAR'RAR  operation;  coherent  and  fully 
polarimetric.  System  parameters:  frequency,  33  GHz;  resolution  (SAR),  1  ftx  1  ft;  beamwidth,  2°;  polarization  isolation.  30  dB;  sensitivity  (SAR 
Mode).  S  N  tOdB  for;  (7km  Range),  a  =  -30dBm^. 


II.  POL.4RIMETRIC  CLUTTER  MODEL 

We  use  a  non-Gaussian  "product''  model  to  characterize  the 
polarimetric  return  from  ground  clutter.  Note  that  with  a 
Gaussian  model,  each  pixel  of  clutter  in  a  spatially  homoge¬ 
neous  region  of  an  image  has  the  same  average  polarimetric 
power.  A  number  of  authors  have  stated  that  it  is  more 
realistic  to  assume  that  ground  clutter  and  sea  clutter,  for 
example,  are  spatially  nonhomogeneous.  A  non-Gaussian 
model  consistent  with  this  more  realistic  assumption  has  been 
proposed  l.‘'-7|,  and,  in  faet,  the  Gaussian  model  is  actually  a 
special  case  of  the  non-Gaussian  model. 

We  assume  that  the  radar  measurement  vector  Y  consists 
of  three  complex  elements:  HH.  HV,  and  V'V.  Therefore. 


jHH\ 

lHH,^iHH„. 

wv  U 1 

HV,  +  y  WV;, 

(1) 

’  vv' 

\  YV,  +  jVV,,  I 

where  HH,  and  HH^,.  for  example,  are  the  in-phase  and 
quadrature  components  of  the  complex  HH  measurement.  Y 
is  assumed  to  be  the  prcjduct  of  a  complex  Gaussian  vector  X 
(representing  the  speckle)  and  a  spatially  varying  texture 
variable  That  is. 

Y  =  VgX.  (2) 

The  vector  X  is  assumed  to  be  circular  complex  Gaussian 
with  a  joint  probability  density  function  (PDF)  of  the  form; 

/(X)=  ^-expf-X’X  'X)  (3) 

v  }ZI 

where  i  =  (XX')  is  the  polarization  covariance  matrix.  The 
vector  X  is  zero  mean  {£(X)  =  0).  The  covariance  matrix  that 
we  use  for  clutter  data  takes  the  following  form  (in  a  linear- 
polarization  basis): 


I  0  pVy 
0  c  0 
p*s/y  0  y 


rr„„  =  £(|WWp).e 


and  p  = 


F.{\HV\') 

■  H(\HH\')'^' 
FAHH  VV') 


F^\VV\-') 
F(\HH\-)  ■ 


r£(|//«|')-£(|V-'T|-')r 


The  assumption  that  Wf  is  uncorrelated  with  HH  and  V'V  is 
not  always  true  (especially  for  man-made  targets  or  for  a 
polarimetric  SAR  with  cross-talk  between  channels);  hi'w- 
ever.  we  have  found  this  assumptitm  is  valid  for  natural 
ground  clutter. 

We  model  the  product  multiplier  g  as  a  gamma-distributed 
random  variable.  This  assumption  is  not  universal;  the  log¬ 
normal  and  Weibull  models  are  al.so  widely  used.  But  if  we 
assume  that  the  gamma  distribution  is  reasonable,  the  PDF  of 
the  product  multiplier  ,i;  is  specified  by 


a  V  p  / 


where  the  parameters  g  and  ('  are  related  to  the  mean  and 
variance  of  the  random  variable  g: 

£( g) =  ge 

f'(g’)  =  g  *'(*'  + 1)  (7) 

With  the  assumption  that  g  is  gamma  distributed,  the  PDF  of 
the  resulting  vector  Y  =  VgX  is  the  modihed  Bessel  function, 
or  generalized  ^-distribution  (.‘i]. 

I  'Y\ 

2  n  I 

/(Y)  =  -TT777— ^  .JTir,.,:-  ■  (b) 

rr  g  1  (e)|Sl  (gYi  ’Y)' 

If  we  set  g  =  (1  /e)  so  that  the  mean  of  the  texture  variable 
is  unity,  then  in  the  limit  as  this  model  reduces  to  the 

Gaussian  model. 


308 


Vol.  4.  .106- .rix  (IW2) 


III.  MINIMUM-SPECKLE  SAR  IMAGE  PROCESSING 

In  this  section  we  consider  how  to  process  the  three  complex 
measurements — HJi.  HV.  and  VV  (i.e.,  the  vector  Y) — into 
pixel  intensity  in  a  way  that  minimizes  speckle.  The  clutter 
product  model  is  used  to  derive  the  optimal  method  of 
polarimetric  speckle  reduction,  which  can  he  interpreted  as  a 
PWF.  Then  the  amount  of  speckle  reduction  that  can  be 
achieved  by  using  the  PWF  is  determined  theoretically. 

The  measure  of  speckle  we  use  is  the  ratio  of  the  standard 
deviation  of  the  image  pixel  intensities  to  the  mean  of  the 
intensities  (.s  ni): 

1  ,9) 

iti  mean(  y) 

where  the  random  variable  y  denotes  pixel  intensity  .  Given 
the  measurements  HH.  HV.  and  VT',  we  wish  to  construct  an 
image  from  the  quadratic 

y  =  YAY  =  gXAX  (HI) 

where  A  is  a  weighting  matrix  that  is  assumed  to  be  Hermitian 
symmetric  and  positive  definite,  thus  keeping  y  positive.  To 
find  the  optimal  weighting  matrix  A*  (i.e..  the  one  that  results 
in  an  image  whose  pixel  intensities  have  the  minimum  po.s.siblc 
.s'/m).  we  use  the  following  results  from  Ref.  8: 

f:{X’AX)  =  tr(5;  A)  =  X  A,  (11) 


(X  A.  )  -  1  (1(.) 

Using  a  L  agrange  multiplier,  ji.  we  minimize  the  uncon¬ 
strained  functional 


VAR(X  AX)  =  tr(i  ' A)' =  Z  a;  (12) 

z  / 

where  E  is  the  expected  value,  tr  is  the  trace.  VAR  is  the 
variance,  and  A,.  A.,  and  A,  are  the  eigenvalues  of  the  matrix 
i  ■  A.  Combining  Eqs  (7)  through  (12)  yields 

( cL  V  -  f  1  =  f  VAR(X'AX)  ^  VAR(g)  ] 

“I  Eqv)  J  U’(y;)  "e-(X’AX)  E-(g)'l 

(13) 


»  I 

Note  in  Eq.  (13)  that  r  is  a  constant.  Therefore,  minimizing 
sIm  is  equivalent  to  minimizing 


(14) 


If  A.  such  that  A'  =  (4,.  A,.  A,|  is  a  minimizing  solution  of 
(14).  then  so  is  aA.  where  a  is  a  real  scalar.  Therefore,  we 
can  minimize  (14)  by  minimizing  its  numerator 


(I-*)) 


with  the  following  constraint  on  its  denominator 


The  solution  derived  above  is  equivalent  to  applying  a  w  hiten¬ 
ing  filter  to  the  polarimetric  vector  Y  before  forming  the 
image  (sec  Fig.  4).  In  the  whitening  process,  the  vector  Y  is 
passed  through  the  filter  5  '  ’  to  obtain 

W-S  '  'Y  =  vgi  '  'X.  (22) 

The  elements  of  VV  are  complex  random  variables  with  equal 
expected  power.  The  covariance  of  W  is  a  sealed  identity 
matrix;  thus  W  is  said  to  be  white.  As  shown  in  Fig.  4.  the 
optimal  solution  to  the  pi'larimetric  speckle-reduction  prob¬ 
lem  is  simply  to  noneoherently  sum  the  powers  in  the  ele¬ 
ments  of  W: 

y  =  \V'W'.  (2.A) 

hence  the  name  polarimetric  whitening  hlter  (PWF). 

The  process  shown  in  Fig.  4  can  also  be  interpreted  as  a 
change  of  polarimetric  basis  from  a  linear  polarization  basis 
(////.  HV.  VV)  to  a  new  ba.si.s  given  by 


In  this  new  basis,  the  three  polarimetric  channels  are  uncorre¬ 
lated  and  have  equal  expected  power.  Thus  the  optimal  way 
to  reduce  speckle  polarimetrically  is  to  sum  the  powers  non- 
coherently  in  these  three  polarimetric  channels. 

We  have  shown  that  the  PWF  algorithm  processes  the 
polarimetric  measurement  vector  Y  in  a  way  that  minimizes 
SAR  image  speckle.  Furthermore,  the  PWF  is  the  maximum- 
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POLARIMETRIC  WHITENING  FILTER 


•  SIMPLE  INTERPRETATION 


WHITENING" 

FILTER 


=.  W  = 


HH 

HV 


vt 

VV  -  p*  yyHH 
VT(l-lpl’) 


UNCORRELATEO  IMAGES 


•  PWF  IMAGE 


y  =  IHHl’ 


IHVl’ 

2 

vv-p*  ^,/7hh 

Iv^l  " 

Vr(i-lPl’) 

Figure  4.  Minimum-speckle  image  processing  Y  is  a  complex 
vector  containing  the  three  linear-polarization  measurements  Using 
the  whitening  filter  gives  a  new  polarization  basis  W  In  this  basis  the 
three  elements  are  uncorrelated  and  have  equal  expected  power 
The  PWF  image,  y,  is  the  noncoherent  sum  of  the  uncorrelated 
Images. 


likelihood  estimate  (MLE)  of  the  spatial  multiplier  e  ^he 
MLE  can  also  be  shown  to  be  an  unbiased,  minimum-variance 
estimator  of  (i.e..  it  achieves  the  Cramer-Rao  lower  bound) 
(see  Appendix  1). 

Next  we  theoretically  determine  the  amount  of  speckle 
reduction  that  cun  be  achieved  by  using  the  PWF.  Although 
the  PWF  solution  is  independent  of  the  PDF  of  the  spatial 
multiplier  g  in  the  product  model,  the  resulting  s  'm  ratio  after 
speckle  reduction  does  depend  upon  /,  (  g).  Thus,  the  s  m  for 
the  PWF  is 


,  .V  I  _  1  /  4  ,  11 

'  m  ^I’wi-  '  V  '  ^  a  ^ 


(25) 


and  the  .s'm  for  a  single-polarimetric-channel  HH  image  is 


HU 


(2b) 


The  i'  parameter  of  the  gamma  multiplier  appears  in  Eqs. 
(25)  and  (26)  because  the  v'm  includes  fluctuations  in  the 
texture  variable  g.  For  an  ideal  speckle-free  image,  in  fact, 
fluctuations  in  the  terrain  reflectivity  across  the  image  are  still 
present,  thus,  its  sim  is  given  by 


ideal 


(27) 


The  i>  parameter  is  closely  related  to  the  log  standard 
deviation  (m,  in  dB)  of  the  texture  component  s>(  the  clutter. 
This  relation  is  derived  in  Appendix  2.  and  values  of  the  v 
parameter  for  clutter  regions  that  have  <r,  of  1  dB  to  3dB  arc 
tabulated.  By  using  Eqs.  (25)  and  (26)  and  the  results  given  in 
Appendix  2,  we  can  calculate  the  reduction  in  the  standard- 
deviation-to-mean  ratio  achieved  with  the  PWF  (relative  to  a 
single-polarimetric-channel  image).  For  clutter  with  a  spatial 
log  standard  deviation  of  1  dB.  the  .s/m  ratio  of  single-channel 
data  is  1.66  times  larger  than  the  s/m  ratio  of  PWF  data.  For 


clutter  with  a  spatial  log  standard  deviation  of  .’dB.  the  s  in 
ratio  of  single-channel  data  is  1 .45  times  larger  than  that  ol 
PWF  data.  In  the  next  section,  these  theoretical  predictions 
will  he  compared  with  measurements  made  on  actual  data 

IV.  SPECKLE  REDUCTION  RESULTS 

In  the  preceding  section,  we  determined  that  the  optimum 
polarimetric  processing  for  speckle  reduction  is  the  PWF.  and 
we  derived  formulas  to  calculate  the  amount  of  speckle 
reduction  achievable  by  using  the  PWF  In  this  section  we 
show  typical  imagery  gathered  by  the  L  incoln  Laboratory 
SAR  and  visually  compare  single-polarimetric-channel  ima¬ 
gery  with  PWF-processed  imagery  ;  then  we  use  actual  clutter 
data  to  calculate  typical  polarization  covariances  of  trees, 
grass,  mixe'd  scrub,  and  shadows;  finally,  we  apply  PWF 
processing  to  these  clutter  data  and  compare  the  actual 
amount  of  speckle  reduction  with  theoretical  predictions 
based  upon  previously  derived  formulas  [Eqs.  (25)  and  (2b)|. 

A.  Typical  SAR  Imagery.  Figure  5  shows  another  S.AR 
image  of  the  golf  course — this  image  was  constructed  by  lirst 
reducing  speckle  polarimetrically  (using  PW'F  processing  while 
preserving  the  1  ft  x  )  ft  resolution)  and  then  using  noncoher¬ 
ent  spatial  averaging  of  the  I  ft  x  1  ft  PW'F  pixel  intensities 
into  effective  I  m  x  )  ni  resolution  pixels.  The  urea  of  the 
image  shown  iti  Fig.  5  is  approximately  5IH)  x  501)  m  Clearly 
visible  in  the  SAR  image  are  the  pond  and  several  sets  of 
trees,  as  well  as  the  putting  green  located  next  to  the  fairway, 
and  a  larger  set  of  trees  located  behiw  the  golf  course.  Figure 
6  shows  an  aerial  photograph  of  the  golf  course  that  was 
imaged  by  the  Lincoln  Lab  SAR  sensor.  Note  that,  although 
the  aerial  photograph  gives  an  excellent  image  of  the  golf 
course  under  conditions  of  good  visibility,  only  the  SAR 
image  would  be  unaffected  by  such  phenomena  as  rain,  cloud 
cover,  or  fog. 

Figure  7  displays  a  SAR  imtige  of  a  highway  overpass  scene 
( 1  ft  X  I  ft  resolution).  Clearly  visible  in  the  SAR  image  are 
the  guardrails  on  each  side  of  the  overpass  and  the  high 
energy  returns  due  to  the  columns  directly  beneath  the  over¬ 
pass.  Figure  S  shows  an  optical  photograph  t>f  the  highwtiv 
overpass  and  we  can  see  these  supporting  structures.  We 
theorize  that  the  road  beneath  the  overpass  and  the  cylindrical 
concrete  pillars  create  a  top-hat  reflector;  multipath  returns 
from  these  lop-hat  reflectors  appear  in  the  SAR  Image  as 
bright  returns,  displaced  in  range  (Fig.  7). 

B.  Comparison  of  HH  and  PWF  Imagery.  Figure  9  displays 
a  SAR  image  (//H  polarization)  of  the  powerlinc  tower  scene 
(1  m  X  1  m  resolution).  Clearly  visible  in  the  upper  and  lower 
portions  of  the  image  are  two  regions  of  trees  separated  by  a 
narrow  .strip  of  coarse  scrub.  Also  visible  in  the  image,  though 
somewhat  faint,  are  four  powerline  towers  positioned  in  the 
scrub  region  (one  pair  of  towers  in  the  upper  left  of  the  image 
and  one  pair  in  the  lower  right). 

Figure  10  shows  the  corresponding  PW'F  processed  image 
of  the  powerline  tower  scene.  Note  that  in  the  PWF  processed 
image,  the  powerline  towers  have  much  greater  intensity  than 
they  had  in  the  single  polarimetric  channel  HH  image  (Fig. 
9).  In  Fig.  11.  we  give  a  graphical  example  of  how  the 
powerline  tower  and  scrub  clutter  distributions  are  changed  by 
polarimetric  pnveessing.  In  particular,  we  show  histograms  of 
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Figure  5.  SAR  image  of  the  golf  course  scene  (1  m  ^  1  m  resolution)  corresponding  to  the  aenal  photograph  shown  in  Figure  6  This  image 
was  formed  by  first  applying  PWF  processing  to  the  1  ft  coherent  data,  then  spoiling  (4x4  noncoherent  averaging)  to  an  effective  1  m 
resolution. 


Figure  6.  Aerial  photograph  of  the  golf  course  in  Figure  5.  This  photograph  was  taken  in  the  spring  of  1989.  Note  that,  unlike  radar  imagery, 
aerial  photography  of  this  quality  can  only  be  taken  under  clear  weather  conditions. 
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Figure  7.  SAR  image  of  a  highway  overpass  (1  ft  x  i  ft  resolution).  This  image  taken  from  a  siripmap  containing  Highway  US  90  in  New  York, 
shows  the  detail  obtained  from  PWF  processing;  note  the  guard  rails  along  both  sides  of  the  overpass 


data  in  the  powerline  tower  region  and  in  the  scrub  region  for 
the  W//  polarization  image  and  for  the  PWF  proccs.sed  image. 
Clearly,  the  histograms  for  the  PWF-processed  data  exhibit 
much  less  intensity  variation  than  do  the  histograms  for  the 
HH  polarized  data.  In  addition,  since  the  PWF  greatly  re¬ 
duces  speckle  in  the  scrub  region,  the  features  of  the  project¬ 


ed  target  shadow  arc  easily  discerned.  Figure  12  shows  a 
1  ft  X  1  ft  resolution  image  of  the  p<iwerline  tower  and  Its 
shadow  projected  onto  the  ground.  Notice  the  intricate 
shadow  structure  in  the  area  below  the  tower  in  this  image, 
and  the  corresponding  physical  structure  of  the  actual  tower 
as  shown  in  Figure  l.f. 


Figure  8.  Optical  photograph  of  one  side  of  the  highway  overpass  (this  photograph  was  taken  in  November  1991 ).  Note  the  placement  of  the 
cylindrical  support  structures  which  form  top-hat  reflectors  with  the  road. 
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Figure  9.  HH  SAR  image  of  the  powetline  towers  in  a  coarse  scrub  region  (1  m  «  1  m  resolution)  Two  pairs  of  powerline  lowers  are  visible  in 
the  scene,  with  one  pair  positioned  at  either  end  of  the  narrow  strip  of  coarse  scrub  running  diagonally  through  the  image 
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Figure  10.  PWF  SAR  image  of  the  powerline  towers  (1  m  x  1  m  resolution).  Note  that  the  towers  have  much  greater  intensity  than  they  have  in 
the  HH  Image  shown  in  Figure  9. 
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PIXEL  INTENSITY  DISTRIBUTIONS 

POWERLINE  TOWER  AREA 
1  ft  X  1  ft  RESOLUTION 


dB 

Figure  11.  Histograms  of  powerline  lower  data  (1ft  “  1ft  res¬ 
olution)  and  scrub  clutter  data  lor  single-channel  HH  image  (solid 
curves)  and  for  PWF  processed  image  (dashed  curves) 


Figure  12.  SAR  image  of  a  single  powerfine  tower  and  its  shadow 
projected  onto  the  ground.  Visible  in  the  lower  portion  of  the  shadow 
is  the  outline  of  the  upper  support  structure  of  the  lower,  including  the 
insulators. 


Figure  13.  Photograph  of  the  powerline  lower  whose  SAR  signa¬ 
ture  and  shadow  were  shown  in  Figure  12 


C.  Polarization-Covariance  Calculations.  I  he  d.ii.i  need  m 
this  experiment  were  eolleeleil  near  Stoekhridee.  m  ihe 
spring  of  IWU  The  scrub  region  loeiilcil  in  Ihe  \ieinit\  iil  Ihe 
powerline  towers  d  ig  1(1)  was  used  ,is  one  of  our  tvpie.il 
clutter  backgrounds— we  calculated  its  polan/alion 
covariance,  i,  .  to  be 
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Note  that  the  polarization  covariance  lor  scrub  clutter  ni.iv  be 
approximated  (quite  accurately)  by  ihe  general  form  de¬ 
scribed  previously  (Eqs.  (4)  and  (.s)). 

Next  we  evaluated  the  polarization  covariances  ol  tree, 
grass,  and  shadow  regions.  For  each  region,  wc  estimated  the 
clutter  polarization  covariance  parameters  rr„„.  e.  y,  and  p  as 
defined  in  Eq.  (S);  these  estimates  are  given  in  Table  I. 
Although  the  quantity  pv  y  in  Table  I  is  complex,  wc  find  for 
natural  ground  clutter  that  the  imaginary  part  is  negligibly 
small. 

The  polarization-covariance  parameters  of  manmade 
targets  are  quite  different  from  those  of  trees,  grass,  and  other 
types  of  clutter.  For  convenience,  we  considered  the  pow- 
crlinc  towers  shown  in  Figs,  d  and  1(1  lo  be  our  hypothetical 


314 


V(.i  4.  .vxi-iis  I  iwyi 


1  abltf  1 . 

klcltct 

l\il.in/alu>n  ci'v.triciiKc 

p.irciinctcf  v 

i>t  (  1  tl 

•  1  ll  1  clultci 

iuhlf  III.  V  fn  Kano'S  thcofv 

ruklikKii 

my  .isui  ymynl 

‘'ui- 

t- 

> 

f  \  > 

I  rccb 

11 

(1  lu 

H  S*t 

St. ruh  II  'Hi 

tl  '14 

Scrub 

U  1'^ 

1  tiS 

M  Nt 

t  ll  .I's's  l  1 

n  ( t  ■■ 

( irjNs 

1 H  J>xr> 

11  )'l 

l  0^ 

0  >  .4 

Shcultivk  ll  's 

1  1  f  >1  ' 

Sh.iJi'w 

n  IMK4 

II  4.^ 

1  IS 

n  4W 

labk-  IV.  Si.iikI.iuI  iIcm.iiumi-  iJHi  .>1  i  I  H  ■  1  in  .Iuim  d.iij 


Uirgcts  ot  mti.ri.si  Uc  i.\tiin.itci.t  thf  pukiri/'aliDii  auariaiicc 
lit  iHir  h\pii|hcliv.il  lariiols  tmm  scroral  huiKircil  hriphi  peaks 
m  the  siemils  ot  eae(i  tower  I  he  \alue  ohlamerl  lor  '  w.is: 


1  S') 


I  (HI  '  pi IHI 

II  211  yll.2s 
II  .s')  •  ;II.II7 


11.211  •  /II  2.S 

I. 4(1  •  /(I. IHI 

II.  Ill  -  /II  1(1 


(I  s')  /(I  117 

mu  ;()  1(1 

(I  ')')  .  /II  IK) 

(2')) 


I  rom  the  al'ose  pola.’';.i;a!i'i)  ein.inanee  m.itns  sse  see  that 
/.(  H\  I  14  I'i  /III,  ).  whieli  implies  that  the  piuserline 
lowers  i;i\e  an  uniisu;ill>  lame  //I -polari/eil  return  Ihis  is 
due  to  the  phssieal  structure  ol  the  towers.  As  shown  in  1  ig 
l.s.  the  tower  Irames  are  reinloreed  with  steel-strut  latliees 
oriented  at  mans  diltereilt  angles  .Also,  trom  the  ahove 
polarization  eovarianee  matrix  we  see  that  //V  is  eorrelated 
with  HU  and  I'V 


D.  PWF  Speckle  Reduction,  l  o  sahdale  the  speekle-redue- 
tion  lormulas  derixed  prexiously.  we  selected  lour  elntter 
regions — shadow,  grass,  seruh.  and  trees  lor  each  region  we 
calculated  standard-dex  lation-to-mcan  ratios  tor  ////.  Il\  .  and 
\  \  data,  and  lor  the  I’WF  data  I  he  s  in  ratios  ol  clutter  data 
gixen  in  I  able  II  show  lower— and  better— numbers  lor  the 
I’Wb  processed  data 

Ihc  singlc-polarimctric-channel  i  ni  ralios  gixcn  in  Table 
II  can  be  used  lo  compute  the  rms  standard-dex  laiion-to-mean 
ratio  tor  each  clutter  region  liq  (2fi)  xxas  used  lo  estimate 
the  .ipproximate  r  lor  each  region,  and  E  q  ( d.x  I  was  use'd  lo 
predict  the  x  m  of  the  PWF  data.  I  able  III  compares  the 
theoretical  predictions  xxith  the  actual  measured  xaliies.  The 
agreement  between  Ihe  theory  and  metisuremenls  is  very 
good — within  .x'V  in  all  eases. 

Although  the  s  m  ratios  given  in  Table  II  clearly  show  that 
the  PWT  reduces  SAK  speckle,  the  more  important  question 
is  whether  the  clutter  log  standard  deviation  has  a  correspond¬ 
ing  decrease,  because  the  log  standard  deviation  directlv 
affects  largel-dcteclion  performance.  Shown  in  Table  IV  are 
Ihe  log  standard  deviations  computed  from  the  1-ft  resolution 
////.  HV.  and  V'V  data,  and  the  1  ft  PWF  data.  The  PWF 
reduces  the  log  standard  deviation  by  tipproximalely  2.1)  to 
2.7  dB  compared  with  single-channel  1  ft  data. 
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E.  PWF  \Aersus  Adaptive  PWF  Comparison.  In  Kci  )  .m 

.ul.ipiive  1‘WI'  .ilgonihm  was  proposed  ,md  .malxzcd  I  his 
algorithm  .idaptixelx  csiiinatcs  polaii/ation  coxan.mccs  ol 
x.irious  regions  ol  clutter  oxer  the  image  .mil  uses  these  s.mie 
estim.ileil  coxan.mccs  to  mmmn/c  llic  speckle  wiihm  e.ieii 
clutter  region.  Since  the  pol.inzation  coxan.mees  ol  glass, 
trees,  and  sh.iilow  regions  were  loiinil  to  he  xeix  smiil.ii 
(  I  able  1).  we  mxestigaled  the  miproxemeiu  m  speckle  lediic 
til..-,  achievable  using  .idaptixe  I’W  1  processing 

We  eoinp.ircil  the  log  si.mtl.ird  dcxi.iiion  ol  the  clultei 
regions  alter  they  were  whitened  with  the  coxaii.iiice  ot  gi.iss. 
versus  the  rc  ult  when  Ibex  were  whitened  with  the  propei 
covariance  (e  g.,  trees  whitened  vvilh  tree  coxari.mecsl,  ,\s  the 
dat.i  in  Table  \  mdic.itc,  Ihc  use  ot  the  ('loper  covaii.mces 
m.idc  virtually  no  change  m  the  log  si.md.ird  dcvi.itioits  Thus, 
because  the  cox ari.inces  ol  Ihc  dilfcreni  regions  are  so  similar, 
lire  extra  compul.ilion  required  tor  ad.q'inc  I’Vk  I  processing 
IS  not  warranted 

F.  Polarimetric  Averaging  VAersus  Spatial  Averaging.  -\s 

mentioned  earlier,  speckle  can  be  reduced  by  noncohereni 
spatial  averaging,  or  spoiling,  of  the  high-rcsolulion  S.AR 
intensity  dal.i:  however,  noncoherent  spatial  aver.igmg  de¬ 
grades  image  resolution.  The  log  standard  deviations  ol  1  It 
and  I  m  PWl-  il.il.i  are  clearlv  superior  to  the  log  standard 
deviation  ot  1  it  .md  1  m  single  polarimclric-channel  data,  as 
is  shown  bv  the  results  presented  m  Table  \  I 

For  grass  regions.  I’WF  data  at  1  ft  resolution  were  mea¬ 
sured  lo  have  a  log  standard  deviation  of  2.11  dli  .At  the  same 
resolution.  Ihc  Hll  data  were  measured  to  have  a  log  standard 
deviation  of  .s  7dB  and  the  H\  data  had  a  log  standard 
deviation  of  s.bdB.  Thus  polarimetric  averaging  improved  Ihe 
results  over  single-channel  data  by  2.7  dB  Noncoherently 
spoiling  the  single-channel  data  to  1  m  resolution  (i  c..  4  >  4 
averaging)  reduced  ir,  by  approximately  .4.7  dB.  I  dB  better 
than  polarimetric  averaging — but  image  resolution  was  sac- 
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X  m  Ratios  ol  1  1  It  ■ 

1  ft  1  clutter  data. 

IfH 

//V 

t'f 

PWl 

Irccs 

1  S') 

l.h9 

1  .4') 

1  H» 

Scrub 

1  4,t 

\.ll 

1  .ts 

(1.94 

Cirass 

1  i; 

[.(»6 

Mb 

Shadow 

(I.W 

0  W 

1.02 

0.69 

Table  V.  St.indard  deviations  |dH|  of  PWF  vs  adaptive  PWt 
data 
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Table  VT.  Standard  deviations  (dB)  of  1-ft-  and  1-ni-rcsolution  data. 
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riticcd  for  the  improvement.  Spoiling  the  PWF  data  to  l-m 
resolution  yielded  the  elutter  log  standard  deviation  of  l.if  dB. 
an  0.7  dB  improvement  over  the  l-m  HH  data,  and  an  O..S-dB 
improvement  over  the  1  m  HV  data.  Similar  results  were 
obtained  for  tree  elutter  and  for  shadows. 

SUMMARY 

We  investigated  the  polarimetric  speckle  reduction  that  is 
achieved  by  using  the  polarimetric  whitening  filter  (PWF); 
this  processing  method  reduces  SAR  image  speckle  without 
degrading  the  spatial  resolution  of  the  image.  Results  ob¬ 
tained  with  actual  SAR  data  show  that  the  PWF  reduced 
speckle — and  it  also  significantly  reduced  clutter  log  standard 
deviation.  At  1  ft  x  1  ft  resolution,  the  log  standard  deviation 
of  clutter  was  reduced  by  approximately  2.5  dB  relative  to  a 
single  polarimetric  channel.  Such  an  improvement  in  log 
standard  deviation  has  been  shown  to  improve  target  detec¬ 
tion  performance  |H,  9j.  PWF  processing  has  also  been  shown 
to  significantly  improve  the  performance  of  clutter  segmenta¬ 
tion  [1()|  and  texture  discrimination  algorithms  (11|. 

Other  researchers  have  applied  PWF.  adaptive  PWF.  and 
other  polarimetric  prtKessing  methods  to  multifrequency 
polarimetric  SAR  data,  A  brief  discussion  of  some  of  this 
important  research  is  given  in  the  following  paragraphs. 

Toma  et  al.  ll.t)  developed  a  method  for  reducing  speckle 
in  SAR  intensity  imagery  using  multifrequency,  fully 
polarimetric  SAR  data:  their  method,  which  is  a  straight¬ 
forward  extension  of  the  PWF.  was  applied  to  4  m  x  4  m 
resolution  data  gathered  by  the  Jet  Propulsion  Laboratory 
(JPL)  airborne  SAR  (14|.  The  three  complex  polarimetric 
measurements  (HH.  HV.  and  VV)  simultaneously  gathered  at 
F-.  L-.  and  C-  bands  were  combined  into  a  d-dimensional 
nicasurement  vector; 

\  =  (HH,,.  HV,,.  VV,..  HH,  .  HV,  .  VV,  .  HH,  .  HV,  .  FT,  ) 

<.W) 

Then  a  9  x  9  whitening  filter  was  applied  to  the  complex  data 
to  obtain  nine  uncorrelated  intensity  images.  Noncoherent 
averaging  of  these  nine  images  produced  a  SAR  intensity 
image  having  an  amount  of  speckle  reduction  equivalent  to 
that  obtained  by  averaging  nine  independent  samples. 

Lee  et  al.  |15J  also  developed  a  method  for  reducing 
speckle  in  SAR  intensity  imagery  using  multifrequency,  fully 
polarimetric  SAR  data;  they  also  applied  their  method  to  JPL 
SAR  data.  The  method  does  not  use  the  complex  HH.  HV. 
and  VV  data,  but  instead  u.ses  the  (xtlarimetric  intensity  data 
(\HH\\  |WVl\  and  |FF|’)  to  achieve  speckle  reduction.  For 
each  single-frequency-band  (P-.  L-.  and  C-).  a  weighted 
linear  combination  of  the  form 

y  =  \HH\' +  k,\HV\- +  k,\VV\'  (.11) 

was  used  to  construct  a  reduced-speckle  SAR  intensity  image. 


y.  I'he  coefficients  k ,  and  k ,  were  selected  to  mmimi/e  a 
mean-square-error  cost  function  (see  Appendix  .A  of  Ref  15 
for  details). 

Since  Eiq.  (.11)  is  a  weighted  average  of  three  correlated 
images,  one  would  not  expect  to  achieve  as  much  speckle 
reduction  ;i,s  that  obtained  by  averaging  three  mdependem 
samples.  I'he  P-.  and  ('-band  intensity  images  produced 
using  Eq.  (11)  were  combined  using  a  weighted  average 
similar  to  Eq.  (11).  This  final  SAR  intensity  image  was 
observed  to  have  an  amount  of  speckle  reduction  equivalent 
to  that  obtained  by  averaging  six  independent  samples.  Recall 
that  the  approach  of  Toma  et  al.  |ll)  provided  an  amount  of 
speckle  rcducin  .;  equivalent  to  that  obtained  by  averaging 
nine  independent  samples. 

The  reason  for  the  difference  in  speckle  reduction  is  that 
Lee,  et  al.  used  HH  and  VT'  data,  whieh  are  correlated,  T  oma 
et  al.  used  a  whitening  filter  to  produce  nine  uncorrelated 
(independent)  samples. 

In  our  studies,  adaptive  PWF  processing  provided  a  margi¬ 
nal  improvement  in  speckle  reduction  (Table  V  shows  these 
results).  This  is  in  agreement  with  the  findings  of  Lee  et  al 
and  Toma  et  al.  However,  adaptive  processing  has  been  found 
to  blur  imagery  ,  especially  near  sharp  contrast  edges  |S) 
(nonadaptive  PWF  processing  produces  very  clear  imagery .  as 
is  shown  in  this  article.  Also,  adaptive  PWF  processing  has 
been  shown  to  produce  degraded  target  detection  perform¬ 
ance  compared  with  nonadaptive  PWF  processing  |8|. 

APPENDIX  1 :  THE  PWF  AS  AN  ESTIMATOR 

In  this  appendix,  we  show  that  for  a  given  polarimetric 
measurement  veetor  V,  the  PWF  is  the  maximum-likelihood 
estimate  (MLE)  of  the  clutter  tevture  parameter  g.  We  also 
show  that  the  PWF  is  an  unbiased,  minimum-variance  es¬ 
timator  (i.e..  it  achieves  the  t'ramer-kao  lower  bound)  of  the 
texture  parameter  g. 

A.  Maximum-Likelihood  Estimate.  In  MLE.  the  parameter 
g  is  treated  as  if  it  were  deterministic  (nonrandom)  but 
unknown.  We  seek  the  value  of  parameter  g  that  makes  the 
observed  vector  Y  most  likely  .  If  we  denote  the  MLE  of  g  for 
a  given  Y  by  g\,| ,  .  then  g\,| ,  is  implicitly  defined  by 

=0  (.12) 

where  pfYlg)  is  the  eonditional  probability  density  function 
(PDF)  of  the  vector  Y  given  g.  This  conditional  PDF  is  easy 
to  evaluate  since,  given  g,  the  veetor  Y  is  complex  Gaussian 
with  mean  0  and  covariance  gi.  Thus 

/’(''!«)=  ^rnL-rexpl-Y’i  'Y/g)  (.11) 

lilS 

where  p  =  the  number  of  polarizations  (HH.  HV.  and  FV'; 
thus  p  =  1  in  our  case).  To  obtain  the  MLE.  w'c  need  to  find 
the  value  of  g  for  which  p(Y|  g)  is  maximum.  Equivalently, 
we  can  find  the  value  of  g  for  which  log  p(Y  |  g)  is  maximum. 

log|p(Y  I  g)l  =  -log(7r''|S|)  -p  log  g  -  Y’S  'V/g  . 

(-14) 
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Vo!.  4.  ( IW2) 


Differentiating  with  respeet  to  g  yields 


VAR(g,„  ,  )  -VAR(  gX  i.  '\ 


(441 


^loglp(Ylg)|=-^  +  '^^’=0 


which  is  precisely  the  PWF  solution. 

Cramer-Rao  Bound.  Next,  we  investigate  the  estimate  , 
and  determine  (1)  that  the  estimate  is  unbiased,  and  (2)  that 
the  estimate  achieves  the  Cramer-Rao  lower  bound  on  var¬ 
iance  of  the  error  (thus  ,  is  an  efficient  estimate  of  the 
texture  parametc  g). 

First  we  demonstrate  that  g^,, ,  is  unbiased.  From  Eq.  (35) 
we  have 


-  \AR(Xl.  'X), 

P' 

Evaluating  the  above,  using  Eq.  (  12)  yields  the  result 

VAR(g,„,)=^.  (45) 

APPENDIX  2:  RELATING  <t^  AND  v 

The  log  standard  deviation  of  the  texture  variable  g  (denoted 
O',  )  is  delined  to  be 

o,  =  \/VAR(l()log,„  g)  .  (dB)  (46) 


g,„,  Y'v  'Y=  i  i  'X.  (36) 

Taking  the  expectation  with  respect  to  X  gives 

'X).  (37) 

But  Eq.  (11)  gives 

£(x'x  'x)  =  frij;  'f:(xx')|  (3k) 

=  fr(i  'l)  =  p. 

Thus  we  have  verified  that  the  MLE  estimate,  g.,,, ,  ,  is 
unbiased.  That  is. 

Next  we  verify  that  g^,, ,  satisfies  the  Cramer-Rao  lower 
bound,  which  states  that  the  variance  of  any  unbiased  esti¬ 
mate.  g  must  satisfy  the  inequality 

VAR(  g)  ’ - ^ - — -  .  (4(1) 

f:j^iogip(Y|g)|} 

C il’itting  the  details  of  the  derivation,  one  may  easily  show 
that 


The  relation  between  o,  and  the  shape  parameter  r  of  the 
gamma  PDF  can  be  described  as 

£'(log  g)  -  [  (log  g)g'  '  exp(  -  I )  </g  .  (47) 

From  Ref  12  |Eq.  (4.352.1)]  we  have 

Edog  g)  =  vEdd log  g  (48) 

(4d) 

Ref.  12  |Eq.  (4.358.2)]  also  gives 

/'.'Klog  g)-']  =  l'P(i')  +  log  g]‘  -t-  f(2.  !■)  .  (.50) 

Thus 

VAR(logg)  =  t(2.  !■)  (51) 

and 

VARdOlog,,,  g)  =VAr(  )  =  (4,34)T(2.  r)  . 

(-“^2) 


f3loglp(Y|,g)]=^  -  .  (41) 

Evaluating  the  expectation  of  the  above,  again  using  Eq.  (11), 
yields  the  result 

.  (42) 


where  from  Ref.  12  JEq.  (9.521.1)] 


f(2,  id= 


I 


„  „  (i’  +  n)' 


(53) 


Table  VII  lists  the  value  of  the  gamma  parameter  as  a 
function  of  the  clutter  standard  deviation  for  some  typical 
clutter  standard  deviations. 


Thus  the  variance  of  any  unbiased  estimator  of  the  texture 
parameter  g  must  satisfy  the  Cramer-Rao  bound 

VAR(g)a^.  (43) 

Finally,  we  can  verify  that  the  unbiased  estimate  g„,,, 
achieves  the  lower  Ixtund. 


Table  VII.  Oamma  parameter  (e)  versus  <>■,  (dB) 


'6 

1.0 

19.3 

l.-S 

8.9 

2.0 

5.2 

2.5 

.3.5 

3.0 

2.6 
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